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ABSTRACT 

In the formation process of black holes, the density and temperature of mat- 
ter become sufficiently high for quarks and pions to appear. In this study we 
numerically investigate stellar core collapse and black hole formation taking into 
account the equations of state involving quarks and/or pions. In our simula- 
tions, we utilize a code that solves the general relativistic hydrodynamics and 
neutrino transfer equations simultaneously, treating neutrino reactions in detail 
under spherical symmetry. Initial models with three different masses, namely, 40, 
100 and 375M , are adopted. Our results show that quarks and pions shorten the 
duration of neutrino emission if the collapse bounces before black hole formation. 
In addition, pions increase the luminosity and average energy of neutrinos before 
black hole formation. We also find that the hadron-quark phase transition leads 
to an interesting evolution of temperature. Moreover, the neutrino event number 
is evaluated for the currently operating neutrino detector, SuperKamiokande, to 
confirm that it is not only detectable but also affected by the emergence of quarks 
and pions for Galactic events. While there are some issues, such as hyperons, 
beyond the scope of this study, this is the first serious attempt to assess the 
impact of quarks and pions in dynamical simulations of black hole formation and 
will serve as an important foundation for future studies. 
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Introduction 



Massive stars with the main- sequenc e mass M > 10 are known to undergo gravi- 
tational collapse at the end of their lives ( IPoelarends et al.l 120081 ) . In particular, stars with 
M < 25M are thought to end their lives as type II supernovae. It is thought that the accom- 
panying explosion is invoked by the shock wave launched by the core bounce due to the nu- 
clear repulsion, leaving a neutron star. On the other hand, the fate of stars with M > 25M© 
can be observationally split into two branches, namely, a hypernova branc h and a faint 



super nova branch, and they are both thought to form black holes eventually ( jNomoto et al. 



20061 ). It has been proposed that strongly rotating massive stars are constituents of the hy- 
pernova branch while nonrotating and weakly rotating massive stars are constituents of the 
faint-supernova branch. Nonrotating stars more massive than faint-supernova progenitors 
are thought to result in so-called failed supernovae, where the shock produced by the bounce 
cannot propagate outward owi ng to the thick outer layer. In fact, a recently discover ed black 
hole candidate with 24-33M (IPrestwich et al.l 120071 ; ISilverman &: Filippenkol 120081 ) may be 
a remnant of a failed supernova, and a survey involving the monitoring of ^10 6 supergiants 



has b een proposed to investigate the end of lives of such massive stars (IKochanek et al. 



20081 ). It should be noted that the hypothesis described above is a hot subject under active 
discussion. 

So far, many numerical simulations of the gravitational collapse of massive stars have 
been performed in an attempt to elucidat e the mechanism of co re-collapse supernovae (e.g., 
Colgate k Whitelll966l : iTotani et al1ll998l : iMarek fc Jankall2009f ). although definitive results 
are still lacking. N umeric al stud ies on black hole formation by stellar collapse have also begun 
to be carried out. iFryerl ( 119991 ) showed that nonrotating stars with 25M < M < 4OM 
produce a faint-supernova explosion. In this case, a proto-neutron star is formed, which 



is th ought to recollapse to a black hole > 10 s after the bounce (e.g., iBaumgarte et al. 



19961 ). On the other hand, nonrotating stars with M > 4OM bec ome failed supernovae 
invol v ing prompt (~1 s a fter the bounce) black hole formation (e.g.. ISumiyoshi et al.l 120071 . 



20081 ; iFischer et al.l 120091 ) . The gravitational collapse of massive stars is accompanied by 



the emission of a large a mount of neutrinos . In the case of SN1 987A, the emitted neutrinos 
were actually detected (IHirata et al.l Il987l ; iBionta et al.l Il987h . Neutrinos are thought to 
play an essential ro le in the explosion mechanism of ordinary core-collapse supernovae (see 
Kotake et al.l 120061 . for a review). Note that, owing to the weakness of their interaction 
with matter, these neutrinos carry detailed information about the dense core that cannot be 
obtained through photons of any frequencies. Incidentally, gravitational waves are another 
candidate to probe the dense core although the detection is challenging. Among various 
types of core-collapse phenomena, the failed supernova is as bright in neutrino emissions as 
ordinary core-collapse supernovae. Since the time evolutions of luminosities and spectra are 
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qualitatively dif ferent from those of supernovae, they can be used to diagnose prompt black 
hole formation ( Sumiyoshi et al.ll2006r iNakazato et al.ll2010 ). 



The collapse of more massive stars with M > 26OM has also been studied assuming 



that they are Population HI (Pop HI) stars, which were t 



re first stars formed in the universe 



Fryer et al.ll2001t INakazato et al. 



12006 



Suwa et ajj 



(e.g. 

be very massive with M > 1OOM (jNakamura fc Umemura 



200 9J). P op III stars are thought to 



200ll ). and they start to collapse 



as a result of the pair-instability during the helium-burning phase. The pair-instability is 
caused by the creation of electron and positron pairs, which consume some of the thermal 
energy to produce the rest masses of electrons and positrons. For the nonrotational case, stars 
with M > 26OM cannot bounce against this collapse, and they form blac k holes directly 
emitt ing a large amount of neutrinos within a shorter (~ 0.1 s) duration (INakazato et al. 



20061 ). On the other hand, stars with 14OM < M < 26OM also undergo pair-instability but 
do not produce black holes. They reverse the collapse by rapid nuclear burning and explod e 
into pieces; these explosions are called pair-instability supernovae (e.g.. IWoosley et al.ll2002l ). 



To investigate the formation process of black holes, the properties of hot and dense 
matter should be considered. In particular, meson condensation, hyperon appearance and 
quark deconfinement are thought to occur at supranuclear density, affecting the equation of 
state (EOS) markedl y. The effects of possib l e phase transitions were once studied for core- 
collapse supernovae (ITakahara fc Satol Il985l : iGentile et al.l 119931 ). T he evolution s of pro to- 
neutron stars including exotic ma tter have been alread y studied by iPons et al.l (j2001aj ) to 
investi gate kaon condensa tion, by iKeil &: Jankal ( 119951 ) to investigate hyperon appearance 
and by IPons et al.l ( 12001bl ) to investigate quark deconfinement. On the oth er hand, prompt 



black hole formation by failed supernovae has been recently s tudied by ISumiyoshi et al. 



( 120091 ) utilizing the EOS including hyperons ( jlshizuka et al.ll2008l ). An EOS including quarks 
and pions was cons tructed and applied t o the collapse of Pop III stars with 1OOM in 
our previous study (INakazato et al.ll2008at ). While it was concluded that quarks and pions 
accelerate the collapse and shorten the duration of neutrino emission, it remains to be seen 
whether these characteristics hold for a wide range of progenitor masses. The purpose of this 
study is to assess the impact of quarks and pions in black hole formation treating general 
relativity and neutrinos for more general cases as the first serious attempt. 

In this paper, we investigate the gravitational collapse and neutrino emission of black 
hole progenitors in a spheric ally symmetric model t aking into account the EOS including 
quarks and pions proposed by INakazato et al.l (j2008al ). Since, as mentioned above, black hole 
progenitors have a broad mass range, we also study the dependence on initial mass. Moreover, 
observational aspects o f emitted neutrinos are addressed following our other previous study 



( INakazato et al.ll2008bl ). This paper is organized as follows. In Section |2l we briefly describe 
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the EOS's, initial models and numerical methods. The main results are reported in Section[3j 
We give the numerical results for the dynamics and neutrino signal of the reference models 
with 40M o in Section 13.14 for the investigation of initial mass dependence in Section 13. 2^ 
and for issues related to neutrino observation in Section 13.31 Finally, Section H] is devoted to 
a summary and discussion. 



2. Setup 



2.1. Equation of State 



In this study, we adopt the EOS formulated by iNakazato et al.l ( l2008af ). which includes 
the hadron-quark phase transi tion for finit e tempe ratures. For the hadronic phase, this EOS 
utilizes a table constructed by IShen et al.l ( 11998aU bl) based on relativistic mean field theory 
with the effects of thermal pions a dded to their table. The MIT bag model of the deconfined 
three-flavor strange quark matter (IChodos et al.lll974j ) is used for the quark phase, and the 
hadron-quark mixed phase is obtained from the Gibbs conditions in the EOS. In this case, 
the substance is composed not only of u quarks but also of d quarks, and there is an essential 
difference fro m the phase transiti on of a single substance such as the liquid-vapor transition 



of H2O (e.g., iGlendenningi Il992[ ). For instance, the pressure of the hadron-quark mixed 



phase is not constant in an isothermal process, as shown in Figure [TJ In the following, we 
refer to the EOS without pions and quarks (the original Shen EOS), the EOS without pions 
and with quarks, the EOS with pions and without quarks, and the EOS with pions and 
quarks as 00, OQ, PO, and PQ, respectively. 

In our hadronic EOS, thermal pions are treated in the minimum model which assumes 
that their effective mass is equal to their rest mass in vacuum. In reality, pions at rest (p = 0) 
are subjected to a repulsive potential in the nucleons and their effective mass becomes larger 
than that in vacuum. In this case, the pion population is suppressed. Thus, our model 
corresponds to an extrem e case where pions are overproduced, provided that the p-wave 
ttN attraction is omitted (jOhnishi et al.ll2009l ). Recently, it has also been p ointed out that 



hype r ons play an important role in black hole formation by stellar collapse (jlshizuka et al. 



20081 ; ISumiyoshi et al.ll2009l ). Hyperons are not included in our hadronic EOS; however, we 
are planning to investigate their effects in a future work. Incidentally, nuclei except for a- 
particles are treated as a single species in the Thomas-Fermi approximation. It is preferable 
to adopt the EOS in nuclear statistical equilibrium (NSE) or, more sophisticatedly, in non- 
NSE abundances determine d from the precedin g quasistatic evolutions, particularly for the 
temperature T < 0.5 MeV jFischer et al.l l2009h . 



-5 - 



In the MIT bag model, free quarks are confined in the "bag", and this bag has a 
positive potential energy per unit volume, B. This parameter is called the bag constant and 
characterizes the model. F or instance, the l arger t he value of B, the higher the transition 



Nakazato et al.l (j2008al ) have shown that the EOS with B 



> 



density and temperature 
250 MeV fm~ 3 (in another unit, B 1 ^ > 209 MeV) is consistent with recent observations 
of compact stars. In the present study, we set B = 250 MeV fm -3 . Note that the end 
point of the transition line (the so-called critical point) has been suggested to reside in 
the temperature range 150 MeV < T c < 200 MeV on the basis of the heavy-ion collision 
experiments and lattice QCD calculations. If this is the case, the quark matter should be 
most stable for all densities at a temperature of >200 MeV and the hadron-quark transition 
should occur at <150 MeV. Our model is consistent with this picture, although it cannot 
reproduce the critical point in principle. We can confirm the validity of our model from the 



phase d iagrams and free energies, which are shown in Figures 5 and 6 of iNakazato et al. 



(I2008al). respectivel y . Ver y recently, dealing with quark matter in the manner described 



above, ISagert et al.l ( 120091 ) performed simulations of the successful supernova explosions of 
10 and 15M progenitors using spherical models. An important difference in their EOS from 
ours is that their bag constant was very small, B = 90 MeV fm" 3 (5 1 / 4 = 162 MeV). Further 
studies are necessary to fix the bag constant. 



2.2. Initial Models 



In this study, we investigate the collapse of stars with various m asses. The adopted 
model s are those of a 40M Q star with solar m etallicity propo s ed by I Woosley fc Weaver 
(119951 ). a Pop III star wit h lOOM^ prop o sed b y iNomoto et al.l (120051 ) and a Pop III star 
with 375M proposed by INakazato et al.l (120061 ) . The former two models were constructed 
from evolutionary calculations, while the latter model was obtained from the equilibrium 
configuration of a gravitationally unstable iron core. The evolution scenario of the model 
with 375M is different from that of the other two models. The 4OM star with solar 
metallicity and the 1OOM Pop III star undergo successive nuclear burning and form an iron 
core at the end of quasistatic evolution. This iron core becomes gravitationally unstable 
owing to the photodisintegration of iron and starts to collapse. These features are similar 
to those in the case of ordinary supernova progenitors. On the other hand, the Pop III star 
with 375M starts to collapse through pair-instability, and an iron core is formed during 
the collapse. There fore, an ad hoc mode l is adopted in this study; however, its validity 
has been veri fied by INakazato et aT. ( 2006) on the basis of the results of recent evolutionary 



calculations ( IFryer et al 



2001 



Ohkubo et al.l 12006( 1 . 
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Note that the collapses of all models ado pted here have alrea dy been examined for the 
EOS model 00, namely, the EOS proposed by lShen et al. (|l998al bh. usine; the same numer- 
ical methods as those utilized in this p a per. The models of st ars w ith 40, 100 and 375M ^ 
were studied by lSumiyoshi et al.l (120071 ) . iNakazato et al. (b()07h . and lNakazato et"aD (j2006l ). 
respectively. Moreover, r esults for the co ll apse o f the 100M© star under EOS models PO, 
OQ and PQ are given in INakazato et al.l ( 12008al ). The collapses of these three progenitor 
models result in black hole formation under EOS model OO, although their dynamical fea- 
tures are different from each other. The 4OM star produces a bounce before black hole 
formation owing to the nuclear force because the central density exceeds the nuclear density. 
The IOOMq star also produces a bounce but owing to the thermal pressure of nucleons. The 
entropy of this model is sufficiently high for nuclei to dissociate into nucleons and a-particles 
at subnuclear density. Then the thermal pressure of nucleons and a-particles produces a 
weak bounce. On the other hand, the 375M star does not produce a bounce and collapses 
to a black hole directly. In this study, we examine whether or not these features are qualita- 
tively and/or quantitatively changed by the effects of quarks and pions through compar isons 
with the results of previous studies ( INakazato et al.ll2006l 120071 ; ISumiyoshi et al.l 120071 ). 



2.3. Numerical Methods 



The general relativistic implicit Lag rangian hydrodynamics code, which simultaneously 
solve s the neutrino Boltzmann equations (lYamada!ll99~7l : lYamada et a l. 1999; ISumiyoshi et al. 
20051 ). is utilized to compute the dynamics of spherical gravitational collapse and neutrino 
transport. This code can solve the evolution of space time as well as the dynamics up to 
black hole formation. Since the event horizon has been proved to always be located out- 
side the apparent horizon, black hole fo rmation can be confirm ed by finding the apparent 
horizon. For the Misner-Sharp metric (IMisner fc Sharp! 1 19641 ). which is spherically sym- 
metric and adopted in our code, the radius of the apparent horizon is simply written as 
r = 2Grh/ c 2 , where c and G are the velocity of light and the gravitational constant, re- 
spectively ( I Van Riperl 119791 ). r and m are the circumference radius and the gravitational 
mass, respectively, and we solve them as functions of time and the baryon mass coordinate. 
Therefore, the appearance of the apparent horizon can be concluded when the coordinate 
at which r = 2Gm/c 2 is satisfied appears in our numerical simulations (see lYamadal Il997l ; 



Nakazato et al.ll2006l for details). 



To obtain the neutrino distribution functions, we solve the Boltzmann equations by a 
finite difference scheme (Sn method) on discretized grid points for the radial Lagrangian 
coordinate, neutrino energy spectrum, and neutrino angular distribution. In our simulation, 



-7- 



we consider four species of neutrino, u e , u e , and 9^, assuming that the distribution function 
of v r (p T ) is equal to that of (i^J. For the collision terms of the Boltzmann equations, we 
calculate the scattering kernels explicitly in terms of the angles and energie s of incoming and 
outgoing neutrinos (see iMezzacappa fc Bruennlll993l ; lYamada et al.lll999l for details). The 
neutrino reactions taken into account are (1) electron-type neutrino absorption on neutrons 
and its inverse, (2) electron-type antineutrino absorption on protons and its inverse, (3) neu- 
trino scattering on nucleons, (4) neutrino scattering on electrons, (5) electron-type neutrino 
absorption on nuclei, (6) neutrino coherent scattering on nuclei, (7) electron-positron pair 
annihilation and creation, (8) plasmon d ecay and creation and ( 9) neutrino bremss trahlung. 
We ad opt the reaction rate for (8) from lBraaten fc Segell (119931) . that for (9) from iMaxwell 



( 119871 ) and that for the others reactions from 



Bruenn 



(jl98a i 



Note that the neutrino treatments described above must be changed for the region where 
the phase transition occurs. This is because our EOS is constructed assuming that electron- 
type neutrinos are in e quilibrium with other p articles in the hadron-quark mixed phase and 
the pure quark phase (INakazato et al.l l2008aj ) . Under this assumption, we set the neutrino 
distribution functions to be Fermi-Dirac functions for all species, and the electron-type 
lepton fraction, Y h is assumed to be conserved for each fluid element, instead of solving the 
transport and reaction of neutrinos to compute the time evolutions of neutrino distribution 
functions and the electron fraction, Y e . Moreover, we neglect the entropy variation resulting 
from the neutrino transport. These modifications are only applied at the mesh points where 
quarks appear. We can justify this simplification in the neutrino treatments because the 
density is sufficiently high for neutrinos to be trapped at the phase transition. In fact, as 
shown later, neutrinos trapped inside the quark core cannot escape because a black hole is 
formed suddenly <1 ms after the phase transition. Note that (3 equilibrium is also used to 
determine the fraction of s-quarks because strangeness is generated only by weak interactions 
such as s <-> u + e~ +u e . If (3 equilibrium is not assumed, the three-flavor quark EOS becomes 
a function not only of density, temperature and Y e (or Yj) but also of strangeness. 

In this study, the numbers of mesh points for the radial Lagrangian coordinate, neu- 
trino energy spectrum and neutrino angular distribution are chosen to coincide with those 
of previous studies given in Section |2"72"1 For instance, for the 4OM model, we use 255 mesh 
points for the radial Lagrangian coordinate and 14 and 6 mesh points for the energy spec- 
trum and angular distribution, respectively, which are the same as those in ISumiyoshi et al. 
(200 7J) . The uncertain ties originating from the resolutions were evaluated to be ~10% by 
Nakazato et al.l (120071 ). Note that rezoning and dezoning of the Lagrangian coordinate are 
performed during the simulations (ISumiyoshi et al.ll2005l ). Since sufficiently high resolution 
is needed in the vicinity of the shock wave, rezoning is performed for the accreting regions. 
On the other hand, dezoning is performed for the inner regions of proto-neutron stars where 
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materials are almost hydrostatic. 



3. Results 
3.1. Collapse of 4OM Star 
In this sect ion, we examine the co llapse of 40M Q star as the reference models. In the 



previous study (jSumiyoshi et al.l 120071 ). the model with EOS 00 was shown to produce a 
bounce owing to the nuclear force before black hole (apparent horizon) formation. Neutrinos 
are emitted mainly during the period from the bounce to black hole formation. In this study, 
we find that these qualitative features are not changed for the models including quarks and/or 
pions. However, quantitative differences appear in, for instance, the time interval between 
the bounce and black hole formation. In the following, we investigate this phase in detail. 

We show the time profiles of the central baryon mass density in Figure [2j The bounce 
owing to the nuclear force corresponds to the spikes at t = 0, which is defined as the time of 
the bounce. At this moment, EOS dependence does not appear because the hadron-quark 
phase transition has not yet occurred and the contribution of pions is still minor. Although 
the bounce produces a shock wave, it does not propagate out of the core and is stalled. At 
the center, a proto-neutron star is formed and gradually contracts owing to the accretion 
of shocked matter. This phase corresponds to the gradual density increase in Figure [2j 
Finally, the core collapses to a black hole and the central density increases rapidly. We 
can see that the time interval from the bounce to black hole formation is reduced owing 
to the contribution of quarks and pions. This is because the EOS becomes softer and the 
maximum mass for the stable configurations of proto-neutron stars decreases as found in 
earlier studies. In fact, the time intervals are 1.049 s, 1.086 s, 1.145 s, and 1.345 s for the 
models with EOS's PQ, OQ, PO, and OO, while the maxi mum masses of the "co ld" neutron 



stars are 1.8M , 1.8M , 2.0M Q , and 2.2M Q , respectively (INakazato et al.ll2008al ). Since the 
mass accretion rate does not differ among the EOS models, a soft EOS leads to a reduction 
of the time interval. Note that we cannot simply quote the maximum masses of cold neutron 
stars because the proto-neutron star is hot and lepton rich. However, our results are roughly 
consistent with this trend. 

While both quarks and pions soften the EOS and promote black hole formation, their ef- 
fects are qualitatively different. Comparing the models with EOS's PO and PQ (or EOS's OO 
and OQ), we can see that quarks have an effect in the very late phase. In other words, the 
transition to the mixed phase triggers the collapse to a black hole. The mass- r adius r elations 



of a neutron star using our EOS's are shown in Figure 7 of INakazato et al.l (j2008al ). From 
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this figure, the maximum mass for the mixed EOS is very close to the mass at which the 
hadron-quark phase transition makes a difference. This description is consistent with the 
fact that black hole formation occurs immediately after the appearance of quarks. On the 
other hand, the effect of pions begins to appear as a gradual increase in the density because 
the thermal pions appear before pion condensation. 

The features of the appearance of quarks and pions can be seen in Figure [3j where the 
profiles of the particle fractions and the baryon mass density of the model with EOS PQ 
are shown for each step. Note that there are nuclei and a-particles in the outer region 
whose profiles are not shown in Figure [3j When the central density is 4 x 10 14 g cm -3 
(t = 272 ms, or 777 ms before black hole formation: upper left panel of Figure [3]), the 
population of thermal pions is small and they do not affect the dynamics. We can also 
confirm this from the comparison of EOS OQ and EOS PQ in Figure [2J When the central 
density is 8 x 10 14 g cm -3 (t = 1022 ms, or 27 ms before black hole formation: upper right 
panel of Figure [3]), pion condensation has already occurred in the central region; however, 
quarks have not appeared yet. When the central density is 2.5 x 10 15 g cm -3 (0.07 ms before 
black hole formation: lower left panel of Figure [3]), quarks in the mixed phase prevail in the 
central region. In this phase, the star is already collapsing dynamically to a black hole, which 
is consistent with the fact that quarks begin to have an effect in the very late phase. At the 
time of black hole formation (lower right panel of Figure [3]), the central density increases to 
1.5 x 10 16 g cm -3 and the pure quark matter resides in the central region. 

We show the time evolutions of the density, temperature, entropy per baryon, electron 
fraction and radial velocity profiles of the model with EOS PQ in Figure HI where the initial 
location of the apparent horizon and the profiles of sound speed with opposite sign at the 
time of black hole formation are also shown. The density rises rapidly after the phase 
transition. In contrast, the temperature profile at the moment of black hole formation has 
a multi-peaked shape. This will be discussed again later. As can be recognized from the 
entropy profiles, the collapse to a black hole is adiabatic. The electron fraction decreases 
during the phase transition owing to the generation of s-quarks, which have negative charge. 
In the pure quark phase, matter is compressed without changing the fractions of not only 
electrons but also other particles owing to the chemical equilibrium. Note that, in this 
regime, quarks and leptons can be regarded as degenerate and relativistic ideal Fermi gases, 
where the number density rii and chemical potential /i« relate as rii oc /i 3 . Comparing the 
radial velocity and sound speed at the time of black hole formation, we can see that the 
infall becomes subsonic for the pure quark phase, while it is supersonic for the hadron-quark 
mixed phase (see also lower right panel of Figure [3]). The sound speed is lower for the 
mixed phase, which corresponds to the gradual increase in the pressure against contraction 
in Figure [lj Therefore, the infall velocity becomes maximum in the region with the mixed 
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phase and does not increase further in the region with the pure quark phase. Note that pure 
quark matter in the MIT bag model is asymptotically close to a relativistic ideal gas at the 
high-density limit. Thus, the sound speed is nearly c/v^3 for the innermost region. 

The evolution of the central density and temperature of the model with EOS PQ is 
plotted with the phase diagram in Figure [5j To draw the phase diagram, the electron-type 
lepton fraction is fixed to Yi = 0.3, which is the same as the value at the center of this 
model. From this figure, we can see that the temperature decreases in the mixed phase 
despite the increase in density. Although th is appears unfamiliar, we can interpret it in 



the context of a phase transition (lMullerlll997l ). Here we assumed that the transition is first 



order, although s ome authors have regarded it as a second-order or crossover transition (e.g., 



Aoki et al.l 120061 ) . In the first-order transition, the release of latent heat occurs. Thus, in 
an isothermal phase transition, the entropy varies. Note that, in our EOS, the low-density 
(hadron) phase has lower entropy, which is opposite to an ordinary liquid-vapor transition 
(e.g., water vapor has a lower density and higher entropy than liquid water in the liquid- vapor 
transition of H 2 0). This means that the entropy increases during the isothermal transition. 
However, in our case, the entropy does not vary because matter in the collapsing star is 
compressed adiabatically. If the onset of the transition point is fixed, the entropy is larger 
for the isothermal transition than for the adiabatic transition (see Figure |6]). Therefore, the 
temperature is lower for the adiabatic transition than for the isothermal transition. This is 
the reason for the temperature decrease during the phase transition. 

The multi-peaked shape of the temperature profile in Figure H] mentioned earlier is due 
to the temperature decrease during the phase transition. For instance, the boundary between 
the hadronic phase and mixed phase is at approximately ~1.4M at the moment of black 
hole formation (lower right panel of Figure [3]). Since the temperature of each Lagrangian 
fluid element (mass coordinate) decreases for the mixed phase, one peak is generated near the 
boundary (see the solid line in Figure H]). On the other hand, a local minimum appears near 
the boundary between the mixed phase and quark phase (~1.1M Q ) because the temperature 
increases again in the pure quark phase. 

We now turn to neutrino emission. In Figure the average energies and luminosities 
of neutrinos are shown as a function of time. Note that and have the same type of 
reactions, the difference in coupling constants is minor and, as already mentioned, v T (z/ r ) is 
assumed to be the same as z/ M ip^). Therefore, we collectively denote these four species as v x . 
The average energy presented here is defined by the rms value. Comparing the models with 
EOS's 00 and OQ or the models with EOS's PO and PQ, we can see that the profiles do 
not significantly differ from each other up to the time of black hole formation for the models 
with quarks. This is because, again, quarks only have an effect at the final moment. As a 
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result, the total energies of emitted neutrinos for the models including quarks are lower than 
those of the models without quarks because of the shorter durations of neutrino emission. 

A similar trend can be seen in the comparison between the models with and without 
pions. However, pions increase the average energies and luminosities of neutrinos gradually. 
Roughly speaking, neutrinos can be regarded as being emitted from the neutrino sphere, 
where the optical depth is 2/3 for neutrinos with a typical energy. The neutrino luminosity 
summed over all speci es is roughly given by the accretion luminosity L^ cc ~ GM V M jR v 



([Thompson et al.ll2003l ). where R u , M, and M v are the radius of the neutrino sphere, the mass 
accretion rate and the mass enclosed by R u , respectively. The average energy of neutrinos 
is approximately proportional to the temperature of the neutrino sphere, T v . Owing to 
appearance of pions, the EOS becomes soft and the inner core contracts as shown in Figure |5J 
Therefore, R v decreases and T v increases. This is why the average energies and luminosities 
of neutrinos are increased by pions. This effect is particularly notable for u Xl because v x 
does not have charged- current reactions. The absence of charged- current reactions makes 
the core optically thinner and the radius of the neutrino sphere smaller for v x . Thus, the 
signal of v x is more sensitive to the difference in the inner region, that is, the appearance of 
pions. 



3.2. Mass Dependences 

In this section, we show the results for the Pop III stars with 1OOM and 375M© 
and compare them with those for the 40M Q star with the solar metallicity given in Sec- 
tio n 13.11 As mentioned alr eady, results for the models with 1OOM have also been reported 



by iNakazato et al.l (j2008al ) . We show the time profiles of the central baryon mass density 
for the models with 100M Q and 375M in Figure El Note that the time is measured from 
the point when the central density exceeds 10 12 g cm -3 for the models with 375M because 
the star collapses to a black hole directly without a bounce. For the models with 100M Q , 
the time is measured from the bounce, similarly to the models with 40M Q . 

The EOS dependence of the models with 100M Q is similar to that for the models with 
40M Q , although there are some differences. First, the central density at the bounce is 
~2 x 10 14 g cm -3 for the 100M o models, while it is ~3.2 x 10 14 g cm -3 for the 4OM models. 
This is because the bounce mechanism is different. As already stated, the 100M o models 
produce a bounce owing not to the nuclear force but to the thermal pressure of nucleons. 
Second, the effect of pions is clearer for the 100M Q models for the following reason. The 
entropy in the central region of a 100M Q star (~3.5A;b per baryon) is higher than that of a 
4OM star (~1&b per baryon), where is the Boltzmann constant, and the temperature 
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is also higher. In Figure [91 we plot the evolution of the central density and temperature for 
the 1OOM model with EOS PQ with the phase diagram for Yj = 0.17, which is the value at 
the center of this model. Upon comparison with Figure one can recognize that a 1OOM 
star has a higher temperature than a 4OM star. Therefore, a 1OOM star has a larger 
population of pions than a 4OM star. Aside from this difference, the fact that quarks and 
pions hasten black hole formation and reduce neutrino emission is unchanged for the models 
with 1OOM . 

For the models with 375M , the EOS dependence is not strong as shown in Figure [SJ 
In this case, the core does not bounce and is already collapsing to a black hole at the time 
of quark and/or pion appearance. Therefore, quarks and pions have a limited effect. In 
Figure [TDl where profiles of the density and radial velocity at the time of apparent horizon 
formation are shown with the locations of the apparent horizon, we can see that quarks 
accelerate the collapse and affect the innermost region inside the apparent horizon. Since, as 
already mentioned, the apparent horizon is always located inside the event horizon, quarks do 
not affect the dynamics and neutrino signals outside the event horizon (black hole). Thus, the 



result s and discussion for the neutrino emissions of Pop III stars reported in iNakazato et al. 
(120061 ) do not require modification. 



3.3. Implications for Neutrino Detection 



When quarks and pions have an impact on black hole formation, the neutrino signals 
detected by terrestrial neutrino detectors will be affected. Among the initial models inves- 
tigated in this study, a difference in the neutrino signal may occur for 4OM and 1OOM 
stars. Our 1OOM star is that of a Pop III star, which is a first-generation star and no 
longer exists survive in the nearby universe. Therefore, in this section, we investigate the 
detectability of neutrinos emitted during black hole formation by our 4OM models for the 
currently operating neutrino detector SuperKamiokande. In the evaluation of neutrin o event 
numbers, we take into account neutrino oscillation following INakazato et al.l (l2008bl ). 



Neutrinos emitted from the stellar core-collapse propagate through the stellar envelope, 
where neutrino flavor conversion occurs by the Mikheyev-Smirnov-Wolfenstein effect. When 
neutrinos pass through the earth before detection, they also undergo flavor conversion inside 
the earth. In this case, the results of flavor conversion depend on the nadir angle of the 
proge nitor. In our a nalyses, we utilize realistic pro files of the progenitor (jWoosley &: Weaver 



19951 ) and the earth (IDziwonski fc Anderson! 1 198 ll ). There are two undetermined parameters 



of neutrino mixing, namely, the mixing angle #13 and the mass hierarchy. For the former, 
only the upper limit is given as sin 2 #13 < 2.0 x 10~ 2 , while the other mixing angles have 



13 



been well measured. Whether the sign of the mass-squared difference Am^ is plus (normal 
mass hierarchy) or minus (inverted mass hierarchy) is also currently uncertain. In this study, 
we investigate the dependence of the neutrino event number on these undetermined mixing 
parameters as well as the nadir angle of the progenitor. The distance from the progenitor, 
R, also determines the event number, which can be scaled simply as oc 1/R 2 . Here we set 
R = 10 kpc, which is a typical length for our Galaxy. 

In Figure [TTJ the time-integrated event numbers of the collapse of the 40M Q models are 
shown for various parameter sets as well as for the EOS models. We can see that the expected 
event number is sufficien tly large for all cases, being similar to that of ordinary supernova 



neutrinos, ~10,000 (e.g.. iTakahashi et al.ll200ll ). The error bars in Figure [TT1 represent the 



upper and lower limits owing to different nadir angles. In the case with t he earth effects, 



the en ergy spectral shape is deformed to a wavelike shape (see Figure 7 of iNakazato et al. 



2008bl . for instance) because the typical length of neutrino oscillation becomes comparable to 



the size of the earth and the neutrino survival probability becomes sensitive to the neutrino 
energy. However, integrating over the neutrino energy, the fluctuation is smoothed out and 
its impact on the total event number is small. For all EOS models, the event number 
decreases for the inverted mass hierarchy with larger sin 2 # 13 because almost all i> e , which 
has the greatest contributions to the event number through the inverse beta decay reaction, 
converts to or v T for this parameter set. 

From Figure [TT] we can see that the emergence of quarks and pions affects the total 
neutrino event number in addition to the duration of neutrino emission for some cases. The 
neutrino event number of the model with EOS OQ is >30% smaller than that of the model 
with EOS 00 for both parameter sets of neutrino oscillation. This difference arises mainly 
from the short duration of the neutrino emission. In fact, as shown in Figured the duration 
for the model with EOS OQ is <20% shorter than that for the model with EOS 00, where 
the luminosity and average energy increase in the late phase. While quarks merely affect the 
duration, pions increase the luminosity and average energy of neutrinos before black hole 
formation. Therefore, the neutrino event number of the model with EOS PO is not greatly 
reduced from that of the model with EOS 00. Moreover, the model with EOS PQ is almost 
identical with the model with EOS OQ. This similarity can be clearly seen in Figure [T21 
which shows the energy spectra of the time-integrated event number. Since the model with 
EOS PQ has a shorter duration but a higher luminosity and average energy of neutrinos 
than the model with EOS OQ, both effects cancel out. 



-14- 



4. Summary and Discussion 

In this study, we have performed a series of black hole-forming core-collapse simula- 
tions involving EOS's with quarks and/or pions. Our EOS's utilize the MIT bag model 
with B = 250 MeV fm~ 3 for the quark phase. We performed numerical computations by 
solving the Boltzmann equation with hydrodynamics under spherical symmetry in general 
relativity to obtain detailed information on the energy spectrum of neutrino emission during 
the evolution. We adopted initial models with three different masses, namely, 40, 100, and 
375M Q , although we mainly reported results for a 4OM model. We found that quarks and 
pions shorten the duration of neutrino emission if the collapse bounces before black hole 
formation. In particular, pions also increase the luminosity and average energy of neutrinos 
before black hole formation. For the Galactic events of ~40M Q stars, the neutrino events 
at the currently operating detector SuperKamiokande are not only detectable but also af- 
fected by the emergence of quarks and pions. The observational features of the model with 
quarks and pions are almost the same as those of the model with quarks but without pions. 
Note that, the hadron-quark phase transition gives rise to the nonmonotonic increase in the 
temperature during the collapse. 

Needless to say, there are some issues beyond the scope of this study. Among other 
issues, there are many other possible EOS's of hot and dense matter. In particular, the 
duration of neutrino emission becomes short and the luminosity and average energy of neu- 
trinos increase also for soft hadronic EOS's. To probe the emergence of q uarks and pions 



for th e Galactic events, more detailed statistical analyses will be needed (INakazato et al. 



20101 ). Hyperons are thought to emerge at the densities of interest, which we have not taken 
into account. While core-collapse simulations invo lving hyperon populat ions have already 
been performed and reported in a separate paper (jSumiyoshi et al. 2009 ). we are planning 



to construct an EOS including quarks, pions and hyperons and apply it to our numerical 
simulations. The consideration of mesons, not only pions but also kaons, requires sophis- 
ticated treatment. The use of a quark EOS with a different value of the bag constant is 
another concern. Moreover, other quark models, such as the Namb u-Jona-Lasinio mode l, 



have been proposed for the modeling of astrophysical phenomena (e.g-. lBlaschke et al.ll2005l ). 
Note that, our result that quarks have an effect in the very late phase may change when 
quark models with lower transition density are adopted. Multidimensional simulations are 
also certainly worth investigating, although their numerical treatment of neutrino transport 
under general relativity is a challenging problem. Uncertainties in the evolutionary calcu- 
lations of progenitor models, such as convect ion and massloss, will affect neutrino sig nals 



through the density profile of the outer layer (jSumiyoshi et al.ll2008l ; iFischer et all 120091 ). It 



should be emphasized, however, that this is the first serious attempt to assess the impact of 
quarks and pions in the dynamical simulations of black hole formation and will serve as an 
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important foundation for future studies. 

We are grateful to Akira Ohnishi and Hideyuki Suzuki for fruitful discussions. In this 
work, numerical computations were partially performed on the supercomputers at Research 
Center for Nuclear Physics (RCNP) in Osaka University, Center for Computational As- 
trophysics (CfCA) in the National Astronomical Observatory of Japan (NAOJ), Yukawa 
Institute for Theoretical Physics (YITP) in Kyoto University, Japan Atomic Energy Agency 
(JAEA) and High Energy Accelerator Research Organization (KEK). This work was par- 
tially supported by Research Fellowship for Young Scientists from the Japan Society for 
Promotion of Science (JSPS) through 18-510 and 21-1189, and Grants-in-Aid for the Sci- 
entific Research from the Ministry of Education, Culture, Sports, Science and Technology 
(MEXT) in Japan through 17540267, 18540291, 18540295, 19104006, 19540252, 20105004, 
21540281 and 22540296. 

REFERENCES 

Aoki, Y., Endrodi, C, Fodor, Z., Katz, S. D. & Szabo, K. K. 2006, Nature, 443, 675 

Baumgarte, T. W., Janka H.-Th., Keil, W., Shapiro, S. L. & Teukolsky, S. A. 1996, ApJ, 
468, 823 

Bionta R. M., et al. 1987, Phys. Rev. Lett., 58, 1494 

Blaschke, D., Fredriksson, S., Grigorian, H., Oztas, A. M. & Sandin, F. 2005, Phys. Rev. D, 
72, 065020 

Braaten, E., & Segel, D. 1993, Phys. Rev. D, 48, 1478 
Bruenn, S. W. 1985, ApJS, 58, 771 

Chodos, A., Jaffe, R. L., Johnson, K., Thorn, C. B. & Weisskopf, V. F. 1974, Phys. Rev. D, 
9, 3471 

Colgate, S. A., & White R. H., 1966, ApJ, 143, 626 

Dziwonski, A. M., & Anderson D. L., 1981, Phys. Earth. Planet. Inter., 25, 297 

Fischer, T., Whitehouse, S. C, Mezzacappa, A., Thielemann, F.-K., & Liebendorfer M. 
2009, A&A, 499, 1 

Fryer, C. L. 1999, ApJ, 522, 413 



- 16 - 

Fryer, C. L., Woosley, S. E., & Heger, A. 2001, ApJ, 550, 372 

Gentile, N. A., Aufderheide, M. B., Mathews, G. J., Swesty, F. D., & Fuller, G. M. 1993, 
ApJ, 414, 701 

Glendenning, N. K. 1992, Phys. Rev. D, 46, 1274 

Hirata K., et al. 1987, Phys. Rev. Lett., 58, 1490 

Ishizuka, C, Ohnishi, A., Tsubakihara, K., Sumiyoshi, K., & Yamada, S., 2008, J. Phys. G: 
Nuci. Part. Phys., 35, 085201 

Keil, W., & Janka H.-Th., 1995, A&A, 296, 145 

Kochanek C. S., Beacom, J. F., Kistler, M. D., Prieto, J. L., Stanek K. Z., Thompson, T. 
A., Yiiksel H. 2008, ApJ, 684, 1336 

Kotake, K., Sato, K., & Takahashi, K. 2006, Rep. Prog. Phys., 69, 971 

Marek, A., & Janka H.-Th., 2009, ApJ, 694, 664 

Maxwell, O. V. 1987, ApJ, 316, 691 

Mezzacappa, A., & Bruenn, S. W. 1993, ApJ, 410, 740 

Misner, C. W., & Sharp D. H., 1964, Phys. Rev., 136, 571 

Miiller, H. 1997, Nucl. Phys. A, 618, 349 

Nakamura, F., & Umemura M., 2001, ApJ, 548, 19 

Nakazato, K., Sumiyoshi, K., Suzuki, H. & Yamada, S. 2008b, Phys. Rev. D, 78, 083014 

Nakazato, K., Sumiyoshi, K., Suzuki, H. & Yamada, S. 2010, Phys. Rev. D, 81, 083009 

Nakazato, K., Sumiyoshi, K., & Yamada, S. 2006, ApJ, 645, 519 

Nakazato, K., Sumiyoshi, K., & Yamada, S. 2007, ApJ, 666, 1140 

Nakazato, K., Sumiyoshi, K., & Yamada, S. 2008a, Phys. Rev. D, 77, 103006 

Nomoto, K., Tominaga, N., Umeda, H., Kobayashi, C, & Maeda, K. 2006, Nucl. Phys., 
A777, 424 



-17- 

Nomoto, K., Tominaga, N., Umeda, H., Maeda, K., Ohkubo, T., Deng, J., k Mazzali, P. A. 
2005, in ASP Conf. Ser., 332, The fate of the Most Massive Stars, ed. R. Humphreys 
k K. Stanek (San Francisco: ASP), 374 

Ohkubo, T., Umeda, H., Maeda, K., Nomoto, K., Tsuruta, S., k Rees, M. J. 2006, ApJ, 
645, 1352 

Ohnishi, A., Jido, D., Sekihara, T. k Tsubakihara, K. 2009, Phys. Rev. C, 80, 038202 

Poelarends, A. J. T., Herwig, F., Langer, N. k Heger, A., 2008, ApJ, 675, 614 

Pons, J. A., Miralles, J. A., Prakash, M. k Lattimer, J. M. 2001b, ApJ, 553, 382 

Pons, J. A., Steiner, A. W., Prakash, M. k Lattimer, J. M. 2001a, Phys. Rev. Lett. 86, 5223 

Prestwich A. H., et al. 2007, ApJ, 669, L21 

Sagert I., Fischer, T., Hempel, M., Pagliara, G., Schaffner-Bielich, J., Mezzacappa, A., 
Thielemann, F.-K., k Liebendorfer M. 2009, Phys. Rev. Lett., 102, 081101 

Shen, H., Toki, H., Oyamatsu, K., k Sumiyoshi, K. 1998a, Nucl. Phys., A637, 435 

Shen, H., Toki, H., Oyamatsu, K., k Sumiyoshi, K. 1998b, Prog. Theor. Phys., 100, 1013 

Silverman, J. M., k Filippenko A. V., 2008, ApJ, 678, L17 

Sumiyoshi, K., Ishizuka, O, Ohnishi, A., Yamada, S., k Suzuki, H., 2009, ApJ, 690, L43 
Sumiyoshi, K., Yamada, S., Suzuki, H., k Chiba, S., 2006, Phys. Rev. Lett., 97, 091101 
Sumiyoshi, K., Yamada, S., Suzuki, H., Shen, H., Chiba, S., k Toki, H. 2005, ApJ, 629, 922 
Sumiyoshi, K., Yamada, S., k Suzuki, H., 2007, ApJ, 667, 382 
Sumiyoshi, K., Yamada, S., k Suzuki, H., 2008, ApJ, 688, 1176 
Suwa, Y., Takiwaki, T., Kotake, K. k Sato, K., 2009, ApJ, 690, 913 
Takahara, M., k Sato, K. 1985, Phys. Lett. B, 156, 17 

Takahashi, K., Watanabe, M., Sato, K. k Totani, T. 2001, Phys. Rev. D, 64, 093004 
Thompson, T. A., Burrows, A., k Pinto, P. A., 2003, ApJ, 592, 434 
Totani, T., Sato, K., Dalhed, H. E., k Wilson, J. R., 1998, ApJ, 496, 216 



Van Riper K. A. 1979, ApJ, 232, 558 

Woosley, S. E., Heger, A., & Weaver, T. A. 2002, Rev. Mod. Phys., 74, 1015 
Woosley, S. E., & Weaver, T. A. 1995, ApJS, 101, 181 
Yamada S. 1997, ApJ, 475, 720 

Yamada, S., Janka H.-Th., & Suzuki, H., 1999, A&A, 344, 533 



This preprint was prepared with the A AS IATgX macros v5.2. 



- 19 - 




Fig. 1. — Pressure as a function of baryon mass density for hadron-quark mixed matter 
(thick solid lines), pure hadronic matter (dashed lines) and pure quark matter (dot-dashed 
lines). The left panel corresponds to the case for matter in neutrino- less /3 equilibrium at zero 
temperature, whereas the right panel corresponds to that with temperature T = 25 MeV 
and electron-type-lepton fraction Y\ = 0.3, where Yi is defined as the sum of the electron 
fraction, Y e , and the electron-type- neutrino fraction, Y Ve . In both the panels, the bag constant 
is chosen as B = 250 MeV fm -3 for the quark matter and thermal pions are included in the 
hadronic matter. 
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Fig. 2. — Time profiles of the central baryon mass density for the collapse of the models with 
4OM . Thin dashed, thick dashed, thin solid and thick solid lines correspond to the results 
for EOS's 00, OQ, PO and PQ, respective ly. The time is measure d from the bounce. Note 
that the result for EOS OO is also given in lSumiyoshi et al.l (120071 ). 
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Fig. 3. — Profiles of the particle fractions (upper plots) and baryon mass density (lower 
plots) of 4OM model with EOS PQ, where Yi = represents the number density of 

particle i, and % represents the baryon number density. The upper left, upper right and 
lower left panels respectively correspond to 777 ms, 27 ms and 0.07 ms before black hole 
(apparent horizon) formation, whereas the lower right panel represents the moment of black 
hole formation. Note that 777 ms before black hole formation corresponds to 272 ms after 
the bounce. 
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Fig. 4. — Time evolutions of the density (upper left), temperature (upper right), entropy 
per baryon (middle left), electron fraction (middle right), and radial velocity (lower left and 
right) profiles for the 4OM model with EOS PQ. Note that plots in the lower-right panel are 
shown as functions of radius while plots in the other panels are functions of the baryon mass 
coordinate. The thick dashed and thick dot-dashed are snapshots 27 ms and 0.07 ms before 
black hole (apparent horizon) formation, respectively, whereas thick solid lines represent the 
moment of black hole formation. Squares show the initial location of the apparent horizon. 
In addition, profiles of the sound speed with opposite sign at the moment of black hole 
formation are shown as thin solid lines in the lower panels. 
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Fig. 5. — Evolution of the central density and temperature of 4OM model with EOS PQ 
(thick solid line) and the phase diagram for Y/ = 0.3. The dashed line represents the 
boundary between hadronic matter and mixed matter, and the dot-dashed line represents 
that between mixed matter and quark matter. The point of inflection at ~3 x 10 14 g cm -3 
is due to the bounce. 
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6. — Schematic picture of the phase transition. 
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Fig. 7. — Average energies (left) and luminosities (right) of neutrinos emitted from 4OM 
models as a function of time after bounce. The panels correspond, from top to bottom, to v e , 
v e and v x (= u^, u T , u^, u T ). Vertical lines represent the end point of the neutrino emission. 
The notation of lines is the same as that in Figure [21 
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Fig. 8. — Same as Figure [2] but for the models with 1OOM (left) and 375 M & (right). Note 
that the time is measured from the point when the central density exceeds 10 12 g cm -3 for 
the models with 375 M , whereas it is measured from the bounce for the models with 1OOM 
and 4OM (Figure |2J. 
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Fig. 9. — Same as Figure [5] but for the evolution of 1OOM model with EOS PQ and the 
phase diagram for Yi = 0.17. 
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Fig. 10. — Profiles of the density (left) and radial velocity (right) for the models with 
375M at the time of apparent horizon formation. The notation of lines is the same as that 
in Figure EJ In the close-up plots in the left panel, the empty triangle, filled triangle, empty 
square and filled square, show the locations of the apparent horizon for the models with 
EOS's 00, OQ, PO and PQ, respectively. 
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Fig. 11. — Time- integrated neutrino event numbers of 4OM models for the normal mass 
hierarchy (left) and the inverted mass hierarchy (right). Error bars represent the upper 
and lower limits owing to different nadir angles. In each panel, the upper left, upper right, 
lower left and lower right plots correspond to the models with EOS's PQ, PO, OQ, and 00, 
respectively. 
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Fig. 12. — Energy spectra for the time-integrated event number of neutrinos in the case 
without the earth effects. The left and right panels correspond to the cases of the normal 
mass hierarchy with sin 2 #13 = 10~ 8 and the inverted mass hierarchy with sin 2 #13 = 10~ 2 , 
respectively. The notation of lines is the same as that in Figure [2j 



